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Abstract 

We have carried out ab-initio calculation and study of structural and elec- 
tronic properties of AgAlM2 (M = S, Se, Te) chalcopyrite semiconductors 
using Density Functional Theory (DFT) based self consistent Tight bind- 
ing Linear Muffin Tin orbital (TB-LMTO) method. Calculated equlibrium 
values of lattice constants, anion displacement parameter (u), tetragonal dis- 
tortion [r] = c/2a) and bond lengths have good agreement with experimental 
values. Our study suggests these semiconductors to be direct band gap semi- 
conductors with band gaps 1.98 eV, 1.59 eV and 1.36 eV respectively. These 
are in good agreement with experimental value within the limitation of lo- 
cal density approximation (LDA). Our explicit study of the effects of anion 
displacement and p-d hybridization show that band gap increases by 9.8%, 
8.2% and 5.1% respectively for AgAlM2 (M = S, Se, Te) due to former effect 
and decreases by 51%, 47% and 42% respectively due to later effect. 

Keywords: A. Chalcopyrite; A. Semiconductors; E. Density Functional 
Theory; E. TB-LMTO 



1. Introduction 

The ternary semiconducting compounds with the formula A^ B^^^C2^ and 
such as CuAlS2, AgAlSe2, CdSiP2, ZnSnAs2 etc have been 
widely investigated because of their potential applications in electro-optic. 
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optoelectronic and nonlinear optical devices. These compounds are promis- 
ing candidates for solar cells [ij , photovoltaic detectors, light emitting diodes 
modulators, filters such as optical light eleminator filters {sl and optical 
frequency conversion applications in solid state based tunable laser systems 
. Among the two chalcopyrite type crystals the most intensive studies have 
been carried out in the series of / — /// — VI2 family compounds where Cu 
is involved as a group I element js], [oj. Veryfew works have appeared in lit- 
erature where the group I element is Ag [71, |8[. No detailed experimental and 
theoretical study of A^B^^^C^^ type semiconductors where group I element 
is Ag and group III element is Al have been carried out. The substitution of 
a group I and a group III elements give an average result similar to that of 
a group II element in II- VI compounds. So in case of AgAlM2 the presence 
of three different kind of atoms provide desired structural, electronic and 
optical properties to satisfy the criteria for suitable electronic devices js], 0] . 
I —III — VI2 type chalcopyrite semiconductors are also important from chem- 
istry and physics point of view in many aspects. Due to presence of group 
IB transition metal the d-electrons participate in chemical bonding. There- 
fore it is important to know how much contribution of d-electrons are there 
and how they influence the electronic properties. Many defect systems of 



this type of chalcopyrites have also been synthesized and studied [10|, lUl, ll2 



Due to defect structure, the system is more porous and so they have rich 
chemistry. Therefore the defect systems can be doped by suitable elements 
to synthesize new materials 13|, [14 . 

Jaffe and Zunger 15| have proposed theory for band gap anomalies for 
^/ j^m^y/ chalcopyrite semiconductors relative to the binary analogous with 
the detailed study of Cu-based chalcopyrites. They have shown significant 
effects of structural distortion on the band gap. Their calculation was based 
on the Potential- Variational-Mixed Basis (PVMB) taking the experimental 
lattice parameters. Similar effects are also expected for Ag-based chalcopy- 
rite but there are no systematic study carried out for these compounds. We 
fill up this gap. We have carried out a detailed band structure calculation us- 
ing localized Tight Binding Linear Muffin- Tin Orbital (TB-LMTO) method 



17|. The ab-initio calculation is based on Density Functional Theory 



(DFT) [18|, |19|. Unlike taking experimental lattice constants we have cal- 
culated total energy to find the equillibrium lattice constants, tetragonal 
distortion {1]) and anion displacement parameter (u). 

Shay and Kasper 20| have shown that for Cu-based chalcopyrite the reduc- 
tion in the band gaps relative to their binary analogous is co-related with the 
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existance of d-bonding. Their calculation was based on the study of spin 



orbit splitting. This is further varified by Jaffe and Zunger [15|. Artus and 
his group [8] have measured spin-orbit interaction at the fundamental gap 
for AgGaSe2- They found that this splitting is smaller than the splitting for 
analogous ideal Zinc blende structure. They argued that this is because of 
p-d hybridization. We have also shown with detailed calculation of partial 
density of states (PDOS) that hybridization of d-orbital of Ag with p-orbital 
of anion have significant contribution to band gap. In TB-LMTO method, 
the basis functions are localized. Therefore, very few basis functions are 
required to represent the highly localized d-orbital of Ag in the systems un- 
der study. Hence the calculation is not only cost effective, it gives also the 
accurate result. 



2. Methodology 

To start with, for our electronic structure calculations we have used the 
well established TB-LMTO method, discussed in detail elsewhere 16|, 111 . 
Electron correlations are taken within LDA of DFT 18, 191. We have used 



the von Barth-Hedin exchange j2ll] with 512 k-points in the irreducible part 
of the Brillouin zone. The basis of the TB-LMTO starts from the minimal set 
of muffin-tin orbitals of a KKR formalism and then linearizes it by expanding 
around a 'nodal' energy point E^^. The wave-function is then expanded in 
this basis : 



La 
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C£ and are TB-LMTO potential parameters and 5'££l (k) is the structure 
matrix. 
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3. Result and discussion 

3.1. Structural properties : 

The tetragonal unit cell of a typical clialcopyrite semiconductor consists 
of two Zinc blende unit cells. There are four Ag atoms, four Al atoms and 




Figure 1: One unit cell of the chalcopyrite lattice 



eight M (M = S, Se, Te) atoms per unit cell. The positions of the various 
atoms in the tetragonal unit cell are : A(Ag) : 0.0 0.0 0.0; B(A1): 0.0 0.0 

0. 5; C(M): u 0.25 0.125, where 'u' is anion displacement parameter. In an 
ideal zinc-blende structure of binary compound like GaAs each anion has 
four similar cations as nearest neighbour. So the bond lengths are equal 
{Rac = Rbc) and the charge distribution is identical around each bond. 
Therefore in binary compound, having zinc blende structure, u is 0.25 and r] 
= 1. We refer this as ideal case. But in chalcopyrite system under study each 
anion has two Ag and two Al cations as nearest neighbours as shown in figure 

1. Due to dissimilar atoms as neighbours the anion acquires an equilibrium 
position closer to one pair of cation than to other. This new position of an- 
ion is called anion displacement which is the main cause for bond alternation 
{Rac 7^ Rbc)- We refer this case as 'non-ideal' case {u ^ 0.25 Sz rj 1). The 
quantity u — ^ = {Rac ~ ^bc)/^'^ measures the extent to which the bond 
will alter in a system. The bond alternation gives rise to structure anomalies 
relative to the ideal zinc blend structure and it has significant effect on band 
gap. 

For self consistent calculation we introduce empty spheres because the pack- 
ing fraction is low due to tetrahedral co-ordination of ions. We ensure proper 
overlap of muffin tin spheres for self consistency and the percentage of overlap 
is found. Table 1 shows the structural parameters including bond lengths of 
the systems under study. Our obtained lattice parameters have good agree- 



ment with experimental results |22| as well as other computational work [15 
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Table 1: Structural parameters including bond lengths of AgAlM2. Og™, Cexp, u^xp are 
experimental parameters [2^ and Uother- other calculated parameters 



Systems 


a 

(A) 


c 

(A) 


(A) 


Cexp 
(A) 


V 


u 


Uexp 


Uother 


Rac 
(A) 


Rbc 
(A) 


AgAlS2 


5.48 


10.90 


5.72" 


10.13" 


0.994 


0.265 


0.290^^ 


0.288" 


2.42 


2.32 


AgAlSe2 


5.78 


11.52 


5.95" 


10.75" 


0.996 


0.263 


0.270" 


0.287" 


2.54 


2.45 


AgAlTei 


6.22 


12.28 


6.29" 


11.83" 


0.987 


0.261 


0.260" 


0.285" 


2.72 


2.64 



" Ref.Ei 



Jaffe and Zunger [15| have calculated 'u' and bond lengths using experimen- 
tal values of the lattice parameters obtained by Hahn 2^ and conservation 
of tetrahedral bonds (CTB) plus 77 = rjtet rule. In our work, we have calcu- 
lated the ground state equilibrium value of lattice parameters, u, 77 and bond 
lengths by minimizing the total energy using first principle procedure. Our 
result shows that u increases from the ideal value u = 0.25 for all the three 
chalcopyrites which is in agreement with the other work of Jaffe and Zunger 



15|. Therefore, we find Raq-m > Rai-m- The x-coordinate of M-atoms is u 
which is a function of lattice constants a and c. Thus the two near-neighbour 
distances A-C and B-C are given by equation 2 and equation 3. 



Rac — « 
Rbc = 0. 



16 
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(2) 
(3) 



3.2. Electronic properties 

(i) AgAlS2 '■ The band structure and total density of states (TDOS) (fig- 
ure 2) show four major subbands of different band widths below fermi level 
Ep within the energy range to -13.84 eV. The first three subbands having 
band widths 2.29 eV, 0.91 eV, 3.17 eV respectively below the valence band 
maxima are separated by very narrow band gaps of 0.22 eV and 0.11 eV. 
The lowest subband having band width 1.24 eV is formed mainly due to 
the contribution of Sulfur 3s states. There is a large band gap of ^ 5.90 
eV between the lowest and third subband. Partial density of states (PDOS) 
for Ag-d and S-p are plotted in figure 3 to study the orbital character and 
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Figure 2: Band structure and TDOS for non-ideal AgAlS2 



nature of hybridization. These figures show that Ag-d and S-p hybridization 

contribute to upper valence band near fermi level. There is a very weak 
contribution of Ag-d states in conduction band. Calculated PDOS for s,p,d 
states of Ag, Al, S show that the contribution to conduction band is due to 
Ag-p, Al-p, S-d and S-p states. 

(ii) AgAlSe2 : Figure 4 for band structure and TDOS show four subbands 
of different band widths below fermi level Ep within the energy range to 
-13.73eV. The lowest subband is formed due to Se-4s states within the en- 
ergy ranges -12.71 to -13.73 eV. There is a large band gap of nearly 6.24 
eV between the lowest and just above subband. The second lowest band 
is mainly formed due to the contribution of Al-s and Se-p states. A very 
narrow band gap of about 0.36 eV separates the second lowest subband from 
the third one. The main contribution to the third lowest subband is due to 
Al-p, Se-p and Ag-d states. The upper most valence band is separated by 
the previous subband by a band gap of about 0.24 eV. Figure 5 for PDOS of 
Ag-d and Se-p depicts that all the Ag character is concentrated at the up- 
per valence band leaving a very negligible amount in conduction band. The 
upper valence band is dominated by Ag-d and Se-p hybrid character. Cal- 
culated PDOS for s,p,d states of Ag, Al, Se show that the conduction band 
has main contribution from Ag-p, Al-p and Se-p and a very less contribution 
from Ag-s and Al-s states. 
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Figure 3: Partial density of states (PDOS) of Ag-d and S-p for non-ideal AgAlS2 



(iii) AgAlTe2 '■ Unlike AgAlS2 and AgAlSe2, figure 6 for electronic band 
structure and TDOS of AgAlTe2 shows three major subbands of different 
band widths below fermi level Ep within the energy range to -12.03 eV. 
The upper most valence band is in the region of to -4.63 eV. It mainly 
consists of Ag-d with some contribution from Te-p states. There is a very 
narrow band gap of about 0.48 eV between the upper most valence band and 
the previous subband. This second subband is in the region -5.11 to -6.24 
eV. It is mainly formed due to the contribution of Al-s and less contribution 
from Te-p states. Whereas the lowest or third subband is formed in the re- 
gion of -11.0 to -12.03 eV due to the contribution of Te-5s states. Figure 7 
for PDOS shows Ag-d character is concentrated in the upper most valence 
band and a very negligible amount of Ag-d character is in conduction band. 
Like in case of AgAlSe2 the main contribution to upper most valence band 
is due to Ag-d and Te-p hybrid orbitals. Calculated PDOS for s,p,d states 
of Ag, Al, Te depict that the main contribution to conduction band is from 
Te-p, Ag-s/p and Al-s/p states. 

In all the cases valence band maximum (VBM) and conduction band mini- 
mum (CBM) are located at center of Brillouin zone denoted as 'G'(r point). 
This indicates that they are all direct band gap compounds having band 
gaps 1.98 eV, 1.59 eV and 1.36 eV respectively. It is known that LDA un- 
derestimates band gap by 30%. If we correct this error, our results are in 
good agreement with experimental results i.e. 3.13 eV, 2.55 eV and 2.27 eV 



respectively for AgAlS2, AgAlSe2 and AgAlTe2 |23|, |24 
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Figure 4: Band structure and TDOS for non-ideal AgAlSe2. 



3.3. Effect of p-d hybridization on electronic properties 

Unlike binary semiconductors like GaAs etc, the ternary semiconductor's 
'A' atoms in B^^^C2^ belong to group I and are transition metals (Cu, 
Ag). So these transition metal's d-orbitals also participate in forming valance 
bands. Therefore we expect Ag-d & anion-p hybridization in the band for- 
mation in ternary semiconductors whereas there is strong s-p hybridization 
in case of binary semiconductors. It is observed in the case of Ag-based 
ternary compounds that p-d hybridization affect the band gap. In order to 
see explicitely the effect of p-d hybridization on the band gap and in general 
on electronic properties, we have calculated band structure and TDOS with- 
out the contribution of the d-orbitals for ideal AgAlM2 systems. Therefor 
we have first freezed the d-electrons and have treated these electrons as core 
electrons. Figures 8,9,10 respectively show the band structure and TDOS 
with Ag-d electron as frozen for AgAlM2. We have summarized the band 
gaps with and without contribution of d-electrons of Ag in table 2. The 
calculated result shows that there is a significant reduction of band gaps due 
to p-d hybridization and it is 51% for AgAlS^, 47% for AgAlSe2 and 42% 
for AgAlTe2. p-d hybridization not only has strong effect on band gap re- 
duction, it also strongly affects the electronic structure in general as seen 
from the comparision of TDOS with (figure 11) and without (figures 8,9,10) 
p-d hybridisation for ideal case. The p-d hybridization can be interpreted 
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Figure 5: PDOS of Ag-d and Se-p for non-ideal AgAlSe-i 
Table 2: Reduction in band gap(eV) due to hybridization for ideal structure. 



Systems 


Hybridization 


Reduction 


With 


Without 


AgAlS2 


1.70 


3.46 


51% 


AgAlSe2 1.42 


2.67 


47% 


AgAlTe2 1.24 


2.15 


42% 



on the basis of simple molecular orbital considerations 15|. The p-orbitals 



that possess the Fis symmetry hybridize with those of the d-orbitals that 
present the same sysmetry. This hybridization forms a lower bonding state 
and an upper antibonding state. The antibonding state that constitutes the 
top of the valence band is predominantly formed by higher energy anion p- 
states. And bonding state is constituted by the lower energy cation d-states. 
Perturbation theory [25| suggests that the two states ri5(p) and ri5(d) will 
repel each other by an amount inversely proportional to the energy difference 
between p and d states. So this rising of the upper most state causes a gap 
reduction. 



3.4- Structural effect on electronic properties 

If we compare the TDOS of ideal AgAlM2 (Figure 11) with that of non- 
ideal TDOS (figures 2,4,6), we observe that the structural distortion like bond 
alternation and tetragonal distortion are important factors which controls the 
band gap. 

The structural contribution to the band gap is mainly controlled by bond 
alternation and has only a small contribution from tetragonal distortion. 
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Similar effect on band gap is seen in CuInSe2 15| . Important thing is to be 



noted that unhke in case of CuInSe2 E15| the bond alternation is responsible 
for the increase in band gap in the case of AgAlM2. Table 3 shows the band 
gap increament in the three cases. The increament is maximum in the case of 
AgAlS2 because anion displacement is maximum (u = 0.265) in this system. 
As u increases, the hybridization between Ag-d & S-p decreases. Therefore 
the occupied valence bands are stabilized by moving down in energy and the 
antibonding conduction bands are destabilized by moving up in energy. Thus 
band gap increases. 

Our study shows that structural distortion not only affects the band gap 
but also the overall electronic structure significantly. A close comparision of 
TDOS for ideal (figure 11) and non-ideal (figures 2,4,6) show distinct differ- 
ences in the structure in DOS. For example a very sharp peak is found at 
energy -3.79 eV for non-ideal AgAlS2 (figure 2) compared to the correspond- 
ing ideal case (figure 11). This is due to individual character of Ag-d and S-p 
states. Figure 2 also shows that second and third bands get separated due 
to distortion in case of AgAlS2 and hence increament of band gap by 0.11 
eV. Similar result is observed in case of AgAlSe2 and AgAlTe2 as shown in 
figure 11. For example a sharp peak at ~ 4.64 eV which appears in ideal 
AgAlTe2 case disappears due to distortion. This is due to weaker p-d hy- 
bridization in non-ideal case. The width of the upper most valence band in 
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Figure 7: PDOS of Ag-d and Te-p for non-ideal AgAlTe2 



Table 3: Increment in band gap(eV) due to structural distortion (without hybridization). 



Systems 


Ideal 


Non-ideal 


Increment 


AgAlS2 


3.46 


3.80 


9.8% 


AgAlSe2 


2.67 


2.89 


8.2% 


AgAlTe2 


2.15 


2.26 


5.1% 



case of AgAlTe2 also decreases by 0.24 eV. There are effects on conduction 
band also. Like in the case of AgAlS2 there is a sharp peak at 5.85 eV which 
disappears due to distortion. The conduction band width also decreases, 0.6 
eV for AgAlS2 and 0.36 cV for AgAlTe2. The effect of distortion on valence 
and conduction bands shows that distortion is also responsible for significant 
change in optical properties of such semiconductors. 

4. Conclusion 

Calculations and study of AgAlM2 (M = S, Se, Te) suggest that these 
compounds are direct band gap semiconductors with band gaps of 1.98 eV, 
1.59 eV and 1.36 eV respectively. Our study further shows that electronic 
properties of these semiconductors significantly depend on structural distor- 
tion and the type of hybridization. The calculation is carried out using DFT 
based TB-LMTO method, we have used LDA for our exchange co-relation 
functional. Taking into account of the underestimation of band gap by LDA, 
our result of band gap and structutal properties agree with experimental val- 
ues. Detail study of TDOS and PDOS shows that p-d hybridization between 
Ag-d and anion-p drastically reduces the band gap. The reduction is 51%, 
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Figure 10: Band structure and TDOS for ideal and without hybridization AgAlTe2- 



47%, 42% respectively for M = S, Se and Te. But band gap increases due to 
anion displacement in contrast to the result obtained by Jaffe and Zunger 15 
in case of CuInSe2- This is due to the relative value of anion displacement 
parameter 'u'. In the case of AgAlM2, u increases whereas for CuInSe2 it 
decreases with respect to the ideal value u = 0.25 corresponding to binary 
zinc blende structure. The gap enhancement due to anion distortion is 9.8%, 
8.2% and 5.1% respectively for AgAlM2 (M = S, Se and Te). TDOS and 
PDOS further show that there is significant effect on the electronic properties 
due to structural distortion and the presence of p-d hybridization. 
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